Bloch Oscillations of an Exciton-polariton Bose-Einstein Condensate 
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We study theoretically Bloch oscillations of half-matter, half-light quasi-particles: exciton- 
polaritons. We propose an original structure for the observation of this phenomenon despite the 
constraints imposed by the relatively short lifetime of the particles and the limitations on the en- 
gineered periodic potential. First, we focus on the linear regime in a perfect lattice where regular 
oscillations are obtained. Second, we take into account a realistic structural disorder known to 
localize non-interacting particles, which is quite dramatic for propagation-related phenomena. In 
the non- linear condensed regime the renormalization of the energy provided by interactions between 
particles allows us to screen efficiently the disorder and to recover oscillations. This effect is useful 
only in a precise range of parameters outside of which the system becomes dynamically unstable. 
For a large chemical potential of the order of the potential's amplitude, a strong Landau-Zener 
tunneling tends to completely delocalize particles. 

PACS numbers: 71.36.+c,71.35.Lk,03.75.Mn 



I. INTRODUCTION 

Cavity exciton-polaritons (polaritons) are mixed 
exciton-photon quasi-particles. A fundamental differ- 
ence with respect to simple photons is that polari- 
tons are interacting both with themselves and with 
other excitations of the surrounding media. As a re- 
sult, they can exhibit Bose-Einstein condensation^ and 
other related fascinating effects, such as superfluidity^, 
vortex formation^, and Josephson oscillations 5 -! 6 -. A 
strong specificity of polariton condensates with respect 
to atomic ones is the possibility of performing reso- 
nant excitation experiments that allow us to achieve a 
very high degree of control on phenomena such as para- 
metric instabilities 7 -^ or bistability^. A second speci- 
ficity lies in the specific spin dynamics and spin struc- 
ture of polaritons^ which has allowed the demonstra- 
tion of interesting phenomena, such as the formation of 
half vorticesii or the demonstration of multistabilityi^. 
However, polaritons typically suffer both from a short 
life time (in the ps range) and from localization induced 
by the presence of structural disorder, which tends to 
limit the observability of propagation effects. Recent 
progress in growth and technology enables the fabrication 
of very high quality samples in which polaritons can live 
up to 50 ps. Such large values have allowed the genera- 
tion of extended condensates with coherence length larger 
than the size of the sample^ 3 -. These characteristics are 
opening a new research field where it will be possible to 
study propagation phenomena of polariton condensates 
in detail. One of these phenomena is Bloch oscillations 
(BOs). A periodic motion of electrons within the first 
Brillouin zone (FBZ) of a semiconductor biased by an 
electric field was predicted to occur about 80 years ago 
by F. Blochii. This phenomenon has been widely theo- 
retically discussed during the twentieth centuryi 5 -. How- 
ever, its experimental observation requires the coherent 
propagation distance to be larger than the oscillation pe- 
riod. This constraint made the observation of electronic 



BOs in bulk semiconductors impossible. It took 60 years 
and required the fabrication of artificial crystals to finally 
observe this phenomeno n 16 ' 17 . 

BOs have been described and observed in three alter- 
native systems showing extended spatial and long tempo- 
ral coherence, namely, coherent light waves propagating 
in photonic crystal a 18 ' 19 , ultracold atoms^ -, and atomic 
Bose-Einstein condensates (BECs) in optical lattice d 21 ' 22 . 
In the linear regime, the electronic, photonic, and atomic 
systems show very similar single-particle behavior. How- 
ever, they differ strongly in their specific non-linear re- 
sponse. In electronic systems, electron-electron interac- 
tion induces some dephasing, leading to the destruction 
of coherent BOs. On the other hand, photons are non- 
interacting particles, and of course, if they propagate in 
a linear media, no peculiar non-linear behavior is ex- 
pected. The situation is radically different for atomic 
BECs in optical lattices where the transport properties 
are strongly affected by the non-linear effects linked with 
the density of condensed particles. The latter topic is still 
widely under discussion, as rich physics emerges out of 
ffc22r2£. Indeed, while an optical lattice constitutes a per- 
fect periodic optical crystal that can lead to long living 
BOs of dilute (ultra)cold atoms^ 6 -, a BEC can experi- 
ence decoherence^ 7 - and dynamical instabilities^^ due 
to non negligible interactions. Moreover, the account 
of disorder leads to phenomena such as Anderson-like 
localizatio n 30 ' 31 or the breakdown of superfluidit y 32 ' 33 . 
The non-uniformity of the potential tends to prevent 
wave function's spreading, while interactions drive the 
system toward superfluidity This interplay has serious 
consequences for the phase diagram of the syste m 34 ' 35 
and most importantly, the disorder induces extra dephas- 
ing which is revealed by a damping of the matter wave's 
BOs 3 ^ 7 .. 

Currently, several works have already been devoted 
to the behavior of polariton condensates in periodic 
lattices 3 ^— , and this topic has started to attract a lot of 
attention. In this paper we theoretically describe BOs of 
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polaritons in a patterned one-dimensional microwirc. In 
both atomic and optical systems, it is possible to design a 
potential that allows short period oscillations (less than 
1 ps for instance in ReffbSl). This cannot short period 
cannot occur in polaritonic systems, where oscillations 
must take place in the strongly coupled part of the lower 
polariton branch. In total, the width of the first Bril- 
louin zone has to be smaller than one fourth of the Rabi 
splitting. This gives values of the order of 1 — 5 meV and 
implies some lower bound for the BO period that should 
be smaller than the polariton lifetime. 
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In the first part of this work, we propose a realistic 
structure in which polariton BOs could be realistically 
observed and show numerical simulations in the linear 
regime. In the second part of the paper, we therefore 
study the effect of structural disorder on BOs of a po- 
lariton condensate. Indeed, an important peculiarity of 
the polariton system is that it is affected by intrinsic 
sample disorder which can be absent in atomic systems 
and which perturbs much less pure photons' propagation. 
As a result, polaritons are often found to be localized 1 -^, 
which evidently can make the observation of propaga- 
tion effects like BOs difficult. In the non-interacting 
regime, we can expect that the polaritons' motion will 
be blurred (dephased) with increasing disorder poten- 
tial. So, strictly speaking, an uncondensed thermal gas 
of polaritons cannot undergo BOs if the disorder is too 
important. Nevertheless, the formation of an interacting 
polariton BEC should permit us to recover an oscillatory 
behavior as a result of a partial screening of the disorder. 
The growth of interaction energy leads the polariton gas 
toward superfluidity 4 ^ at low velocities, where the influ- 
ence of disorder is the strongest. 



In general, the oscillations are found to be damped be- 
cause of both the residual scattering by disorder when 
polaritons are moving at a supersonic velocity, and the 
occurrence of specific parametric instabilities, discussed 
in Refill and well known in the atomic BEC field, when 
polaritons are accelerated to the inflexion point in the 
first Bloch band of the dispersion. The most interest- 
ing result is that one of these effects can compensate the 
other: The interactions protect (at least partially) the 
BOs from the damping linked with disorder. The precise 
characterization of these non-linear scattering processes 
is not in the scope of the present work. However, we 
would like to point out a key specificity of the polariton 
system with respect to the atomic one: the possibility 
to perform resonant excitation experiments (with spatial 
and frequency resolution) that allow us to address specif- 
ically the various non-linear scattering processes that in- 
teracting particles accelerated in a disordered periodic 
potential can undergo. 



FIG. 1. (Color online) (a) The structure: laterally narrowed 
wire shaped microcavity with periodic metallic depositions 
and the resulting total potential (Only the central part where 
oscillations take place is represented for clarity), (b) The 
associated initial dispersion of the particles modified by the 
periodic potential U ex (x) created by the metal. 



II. STRUCTURE AND LINEAR BLOCH 
OSCILLATIONS 

The structure proposed is schematically shown in 
FigfjTa). It is based on a L x = 100 /im long GaAs micro- 
cavity etched in the y-direction in order to realize a wire 
having a lateral size in the fim range. The confinement 
energy provided by the lateral etching is approximately 
E c = h 2 Tr 2 /2m*Ly, where m* ~ 2m p h (at zero detuning) 
is the effective mass of the polariton. The linear po- 
tential ramp needed in order to mimic the electric field 
acting on electrons can be realized by changing the lat- 
eral size of the wire along its main axis x with a square 
root dependency L y (x) ~ Lq/^/x in the region where 
BOs are expected. In principle, any type of accelerating 
potential can be designed thanks to modern lithography 
technics but of course, there are some limitations due 
to the confinement energy dependance. It is also possi- 
ble to make use of the wedge character of microcavities 
along the growth z-directionii. It should be noted that 
this potential ramp will be acting on the photonic part 
of the quasi-particle. The periodic potential required to 
open a gap in the polariton dispersion can be obtained 
either by depositing a metallic pattern along the wire 4 ^, 
by excitation of a surface acoustic wav o 41 ' 45 , by a square- 
wave-like lateral etching, by a chain of coupled Josephson 
junctions or even by means of Tamm plasmons 4 ^. The 
first two techniques allow the creation of a periodic po- 
tential acting on the excitonic part of the polariton, and 
in this work we will focus on the first one. 

To describe the dynamics of the system we use a 
set of one-dimensional time dependent Scrodinger and 
Gross-Pitaevskii equations for the photonic ip p h{x, t) and 
excitonic fields ip ex (x, t), coupled by the light- matter 
interactional £Ir = 14 meV (Rabi splitting of realistic 
samples), taking into account the finite lifetimes r p h = 40 
ps and T ex = 150 ps of photons and excitons, and a quasi- 
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resonant Gaussian photonic (laser) pump P (t): 
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where m p h = 5 x 10~ 5 m e i, m ex = 0Am e i and m e i are the 
photon, exciton, and free electron masses, respectively. 
The nonlinear term with a = &Ef,a 2 B /S, where E^ is the 
exciton binding energy, its Bohr radius, and S is a 
normalization area, accounts for the polariton-polariton 
coherent repulsive interactions^. Finally, U p h(x) = —Fx 
and U ex {x) are the accelerating ramp potential producing 
a constant force F = 0.2 meV/ p,m and the squarewave 
potential of period d = 1.56 pan and amplitude A = 2 
meV , respectively . Here we make the following remark: 
The barrier height is another important limitation to our 
structure. Indeed, contrary to optical potential where it 
can be easily tuned in wide ranges, for polaritons metal- 
lic depositions will tend to absorb photons, which can 
shorten the polariton lifetime even more. Thus, we can 
hardly consider large enough band gaps to allow us to em- 
ploy the tight binding approximation, which justifies the 
numerical considerations below. Figure [UL) shows the 
modified dispersion of exciton-polaritons with gaps both 
in the upper and lower polariton branch (LPB), opened 
by U ex . In the following, we will concentrate on the LPB, 
where the condensation usually takes place. The gap and 
the first band where oscillations are expected have width 
of A g = 0.75 meV and Ai = 1 meV, respectively. 

To begin our analysis and demonstrate the possibility 
of obtaining regular single-particle BOs in such a sys- 
tem, we show in FigJ5fa) the propagation of a 2-ps-long 
and 20-/im-large photonic Gaussian pulse tuned close to 
the energy of the LPB at k = 0, with an amplitude low 
enough to assume a linear regime. Under the action of 
the constant force and for the parameters we use, our sys- 
tem exhibits BOs of amplitude Abo ~ 12 pm and period 
Tb « 25 ps; the latter is close to lifetime of the particles. 
A further increase in F, would enhance the Landau- Zcner 
tunneling (LZT) probability [visible in FigJ2Ja)] given 
by Plzt — exp (— A 2 m*ir/2F) , and induce a significant 
splitting of the wave function at each period of oscilla- 
tions. Figure [21b) shows the time evolution of the wave 
vector k x along the wire with characteristic reflections on 
the FBZ's edges, which has an extension Zb = 2ir/d = 4 
HTn~ , changing k x into — k x . We emphasize that the 
long polariton lifetime obtained in modern structures is 
crucial for the observation of polariton BOs. With the 
realistic limitations on the ramp potential, it is difficult 
to reduce the period of oscillations below 20 ps. There- 
fore, the polariton lifetime should exceed this value, for 
the phenomenon to be observable. 




FIG. 2. (Color online) Propagation of a photonic Gaussian 
pulse in the biased square-wave potential, (a) The density 
probability versus time and space reveals clear BOs with a 
period of 25.3 ps and an amplitude of 11.7 /jm, and b) density 
vs time and momentum shows the characteristic reflections 
on the edges of the FBZ where the particles are localized. In 
both plots a LZT is visible as a fraction of the non-oscillating 
particles that escape at each period. 



III. BLOCH OSCILLATIONS OF THE 
POLARITON BEC 

In the second part of this work, we consider the real 
specificity of the polariton condensate with respect to 
the coherent photonic wave, namely, the role played by 
interactions. This role is actually better revealed by con- 
sidering the impact of a realistic structural disorder on 
polariton BOs. This disorder can have several origins, 
which could be the intrinsic sample imperfections, the 
non-ideal lateral etching, or some natural (artificial) fluc- 
tuations of the wells' depths and widths in the periodic 
potential. We assume that initially, before the potential 
ramp is applied, the resonantly created polariton conden- 
sate is at thermodynamic equilibrium at T = K, and 
is therefore in its lowest energy state. The ground state 
is found by minimizing at fixed density the free energy 
of the coupled exciton-photon system^: 



E 



dx i j={ph,ex} 



(3) 

Next we solve the time-dependent Eqs.(l) and (2), start- 
ing from the initial condition given by the ground state 
found previously, taking into account the action of the ac- 
celerating ramp potential for different chemical potentials 
(i (given by the energy of the ground state). To clarify, 
we first assume an infinite lifetime of the particles. For 
a particular modulated barrier height of 4 meV, to re- 
duce the LZT probability, which decreases the number of 
oscillating particles, we show different accessible regimes 
(described in the next paragraph) . Figures [3] and @] dis- 
play the probability density of the photonic wave function 
versus time and momentum and the average velocity of 
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FIG. 3. (Color online) BOs in the presence of structural dis- 
order, (left panel) Probability density of the photonic wave 
function versus time and momentum, (right panel) Average 
velocity of the condensate versus time, (a) and (b) Non- 
interacting damped case (/i = 0) for Ud = 0.1 meV (solid 
blue curve) and Ud = 0.2 meV (dotted red curve), (c) and 
(d) Interaction-induced revival of the oscillations for fi — 0.3 
meV. 



the condensate as a function of time. The latter is given 
by 

(v(t)) - ih L [ dx {j>(x,t)Vr(x,t) 1 
{V W) ~ 2m* N (i) J dX \ -V* (x, t) (x, t) / W 
o 

where N(i) the total number of particles. The latter is 
time dependent to compensate for absorbing boundary 
conditions or the lifetime of the particles. Indeed, (v(t)) 
is relevant for imaging the global motion of the conden- 
sate, changing its sign (v (t)) > (< 0) when particles 
change propagation direction. 

The ground state of a non-interacting or, equivalently, 
a very dilute condensate in the presence of disorder is 
mainly localized in the lowest well of the sample. As a 
result, while the condensate is put into motion by the 
constant force, it is very sensitive to the disorder; the 
latter induces back scattering, which leads to a dephas- 
ing and thus either to a damping of the oscillations or to 
their total destruction if the inhomogeneity is strong. In 
that latter case the system as a whole is insulating. In 
real systems we can expect the disorder induced by the 
sample imperfections to be weak with respect to the total 
periodic potential's amplitude. In order to introduce im- 
perfections, we modulate the square-wave potential ran- 
domly in amplitude along the wire with a standard de- 
viation Ud = 0.1 meV from the original potential. In 
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FIG. 4. (Color online) (a)-(b) Unstable and (c)-(d) strong 
LZT regimes. For fj, = 0.5 meV and 3 meV, respectively 



FigsEta) and[3jb) we present this non-interacting case, 
which corresponds to a zero chemical potential (N = 1 
particle). The dephasing is visible in FiglSJa) where the 
creation of new harmonics blurs the oscillations, resulting 
in a damping and deformation of the motion in Figj3th). 
The solid blue curve corresponds to the parameters given 
here while the dottcd-rcd curve is for a disorder that is 
twice stronger, showing the total suppression of the os- 
cillations. 

With the increase in the chemical potential, interac- 
tions drive the system toward superfluidity in the sense 
we used in Refl43l. that is, particles are no more affected 
by the presence of an in-plane potential so they can freely 
propagate in space without being scattered, and the in- 
teraction energy is able to efficiently screen the disor- 
der. This phenomenon is called a "dynamical screening 
of disorder" in RefUfl As a result, the damping is signif- 
icantly reduced. Thus, owing to the interacting nature of 
the particles involved, BOs of exciton-polaritons should 
be observable despite the presence of a structural disor- 
der. Such a situation is depicted in FigOJc) and(3Jd) for 
/x = 0.3 meV . 

Nevertheless, interactions are not sufficient to re- 
cover perfect oscillations obtained in linear homogeneous 
regimes, as one can see in Figj3jd), and they also have 
drawbacks: Indeed, while the condensate is accelerated 
up to the FBZ, at a certain point, the velocity exceeds a 
critical value, and the particles enter a supersonic regime 
where they are no longer superfluid and are thus scat- 
tered on the disorder. Scattering results in a dephasing 
of the oscillations while approaching the FBZ's edges. 
Moreover, with a further increase in the density, para- 
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FIG. 5. (Color online) Situation taking into account the finite 
lifetime of the particles for /i = 0.35 meV. (a) Oscillations 
remain visible, but (b) they are damped due to the decreasing 
interaction energy. 



and losses after a stabilization time. Then we suppose 
that the condensate is loaded by some means in the dis- 
ordered tilted periodic structure while the pump is turned 
off. The condensate will then oscillate a few times and 
exponentially die out. So we, of course, expect the inter- 
action energy to vanish with the particles; thus, a con- 
densate prepared with a sufficient chemical potential to 
screen the disorder will suffer from growing localization 
with decreasing density, which will damp the oscillations. 
So here the ratio between lifetime and oscillation period 
is even more crucial. FigureJS] illustrates such a situa- 
tion, where \i = 0.35 meV is taken slightly higher than 
in FigslSjc) and [3jd) to compensate as much as possi- 
ble for the localization for first oscillation periods. While 
the damping remains significant as seen in [SJb) , several 
oscillations should be observable as shown in[SJa). 



metric instability develop o > ' . This kind of process 
is well known in the field of polaritons££ as the lower 
branch of their bare dispersion possesses an inflection 
point; thus above some density threshold, two polaritons 
can scatter at this point toward signal and idler states, 
conserving both momentum and energy. Here the same 
kind of phenomenon can take place, independently of the 
disorder, within the Bloch bands, which also possess in- 
flection points. We show in FigfJIa) and0|b) this kind 
of phenomenon appearing for \i = 0.5 meV] oscillations 
are completely deformed and no more quantifiable. It 
is important to note that for stronger disorders, it is not 
possible at all to find a screened regime [FigJ3Ja)-(b)], be- 
cause the system enters the parametric instability region 
before recovering oscillations. 

Finally, an even stronger density wipes out all the ef- 
fects mentioned above and puts the Landau-Zener tun- 
neling into play. The interactions renormalize the lowest 
Bloch band, closing the gap in the dispersion and in- 
creasing drastically the probability of the tunneling to 
the second band^. As a result, BOs vanish, and parti- 
cles are almost uniformly accelerated. They practically 
no longer feel the lattice. This regime is illustrated in 
FigsfJJc) and BJd) for a very strong, although probably 
not experimentally accessible, chemical potential of fi = 3 
meV. The condensate is accelerated to the edge of the 
sample and disappears on the absorbing boundaries, as 
seen in[Hc). The particles are no longer reflected at the 
FBZ edges. In|4jd), the mean velocity keeps increasing 
until particles reach the boundary, and then it decreases 
because of the remaining trapped particles. 

To conclude, we will now take into account the finite 
lifetime of the particles to move to a realistic situation. 
We assume that the condensate is created by a quasi- 
resonant pumping close to the LPB at k x = that mas- 
sively populates the ground state of our system. The 
chemical potential is defined by a balance between gains 



IV. CONCLUSIONS 

In conclusion, we have proposed a realistic structure in 
which BOs of exciton-polaritons, a phenomenon not pre- 
viously reported, could be observed. We have analyzed 
both linear and non-linear behavior and found a regime 
where oscillations are able to overcome the disorder ex- 
pected in realistic systems as a result of to the interacting 
nature of the particles. At higher chemical potential, dy- 
namical instability switches on and destroys oscillations, 
and at even higher density the massive Landau- Zcner 
tunneling leads to a strong delocalization of the parti- 
cles. 

Regarding a real experimental investigation, we would 
like to point out some phenomena to be envisaged. First 
of all, as we mentioned in Sec. I, metallic deposition will 
tend to reduce the cavity photon's lifetime; thus a lateral 
(along the x-axis) square-wave etching could be consid- 
ered as a serious option. In a real two-dimensional system 
one should expect transverse excitations to modify the 
instability threshold^. Finally, in a configuration where 
Landau-Zener tunneling is significant (for a weak lattice 
or strong potential ramp) the wave function will split 
and the non-oscillating part will be backscattered at the 
wire's edges, which will slightly blur the motion. Nev- 
ertheless, the exceptional progress in growth and tech- 
nology should lead very soon to even higher-quality sam- 
ples, in which the cavity photon lifetime could approach 
fOO ps. In such samples, Bloch-oscillations of exciton- 
polaritons could be observed, perhaps even in the linear 
regime if the structural disorder is sufficiently low. 
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